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\ We use the high energy predictions of a QCD-inspired parameterization 

of all accelerator data on forward proton-proton and antiproton-proton scat- 
tering amplitudes, along with Glauber theory, to predict proton-air cross 
sections at energies near y/s ~ 30 TeV. The parameterization of the proton- 
proton cross section incorporates analyticity and unitarity, and demands that 
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qj , the asymptotic proton is a black disk of soft partons. By comparing with 

the p-air cosmic ray measurements, our analysis results in a constraint on the 
inclusive particle production cross section. 
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Cosmic ray experiments measure the penetration in the atmosphere of particles with en- 
ergies in excess of those accelerated by existing machines — interestingly, their energy range 
covers the energy of the Large Hadron Collider (LHC) and extends beyond it. However, 
extracting proton-proton cross sections from cosmic ray observations is far from straight- 
forward Jl|. By a variety of experimental techniques, cosmic ray experiments map the 
atmospheric depth at which cosmic ray initiated showers develop. The measured shower 
attenuation length (A m ) is not only sensitive to the interaction length of the protons in the 
atmosphere (A p _ a i r ), with 
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but also depends on the rate at which the energy of the primary proton is dissipated into 
electromagnetic shower energy observed in the experiment. The latter effect is parameterized 



in Eq. (0) by the parameter k; m p is the proton mass and o~ p n ® l a i T the inelastic proton-air 
cross section. The value of k depends on the inclusive particle production cross section in 
nucleon and meson interactions on the light nuclear target of the atmosphere and its energy 
dependence. We here ignored the fact that particles in the cosmic ray "beam" may be 
nuclei, not just protons. Experiments allow for this by omitting from their analysis showers 
which dissipate their energy high in the atmosphere, a signature that the initial energy is 
distributed over the constituents of a nucleus. 

The extraction of the pp cross section from the cosmic ray data is a two step process. 
First, one calculates the p-air total cross section from the measured inelastic cross section 
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Next, the Glauber method is used to transform the measured value of & p n * l ail into a 
proton-proton total cross section a pp ; all the necessary steps are calculable in the theory. In 
Eq. (H) the measured cross section for particle production is supplemented with a p _^ T and 
°p-airj the elastic and quasi-elastic cross section, respectively, as calculated by the Glauber 
theory, to obtain the total cross section <Tp_ air . The subsequent relation between cr^f^r an d 
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<jpp involves the slope of the forward scattering amplitude for elastic pp scattering, —^ z , 
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and is shown in Fig.|T], which plots B against cr pp , for 5 curves of different values of a 1 
This summarizes the reduction procedure from cr^ ir to a pp |[J. Also plotted in Fig. [I] is a 
curve of B vs. a pp which will be discussed later. 

A significant drawback of the method is that one needs a model of proton-air interactions 
to complete the loop between the measured attenuation length A m and the cross section 
°p^air> ^- e v the value of k in Eq. (fj). A proposal to minimize the impact of theory is the 
topic of this letter. We have constructed a QCD-inspired parameterization of the forward 
proton-proton and proton-antiproton scattering amplitudes which is analytic, unitary 
and fits all data of a to t, B and p, the ratio of the real-to- imaginary part of the forward 
scattering amplitude; see Fig.|2|. We emphasize that all 3 quantities are simultaneously 
fitted. Using vector meson dominance and the additive quark models, it accommodates a 
wealth of data on photon-proton and photon-photon interactions without the introduction 
of new parameters M. Because the model is both unitary and analytic, it has high energy 
predictions that are essentially theory-independent. In particular, it also simultaneously fits 
(Tp P and B, forcing a relationship between the two. Specifically, the B vs. a pp prediction of 
the model is shown as the dashed curve in Fig. |I[ The dot corresponds to our prediction of 
o~ pp and B at y/s = 30 TeV. It is seen to be slightly below the curve for 490 mb, the lower 
limit of the Fly's Eye measurement, which was made at y/s ~ 30 TeV. The percentage error 
in the prediction of a pp at \fs = 30 TeV is ~ 3.3%, due to the statistical error in the fitting 
parameters (see references 0, 0]). 

In Fig. ^ we have plotted the values of a pp vs. o~ p n ^l ir that are deduced from the intersec- 
tions of the B-a pp curve with the o~ p tlcA aiT curves of Fig.[I|. Figure [5] allows the conversion of 
the measured cr™^ r to a pp . The percentage error in ct^H^t is ~ 1.9 % near cr^^ ir = 450mb, 
due to the error in a pp from the model parameter uncertainties. 
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Our prediction for the total cross section function of energy is confronted with 

all of the accelerator and cosmic ray measurements in Fig.|j. For inclusion in Fig.|], 
we have calculated the cosmic ray values of a pp from the published experimental values of 
°iJ^airj usm g the results of Fig.|3|. We note the systematic underestimate of the cosmic ray 
points, roughly about the level of one standard deviation. 

It is at this point important to recall Eq. ([!]) and consider the fact that the extraction 
of cri^air fr° m the measurement of A m requires a determination of the parameter k. The 
measured depth X max at which a shower reaches maximum development in the atmosphere, 
which is the basis of the cross section measurement in Ref. 0, is a combined measure of 
the depth of the first interaction, which is determined by the inelastic cross section, and of 
the subsequent shower development, which has to be corrected for. The position of X max 
also directly affects the rate of shower attenuation with atmospheric depth which is the 
alternative procedure for extracting cr^f^. 

The model dependent rate of shower development and its fluctuations are the origin 
of the deviation of k from unity in Eq. (p]). Its values range from 1.5 for a model where 
the inclusive cross section exhibits Feynman scaling, to 1.1 for models with large scaling 
violations |L|. The comparison between data and experiment in Fig. ^ is further confused by 
the fact that the AGASA H and Fly's Eye || experiments used different values of k in the 
analysis of their data, i.e., AGASA used k = 1.5 and Fly's Eye used k = 1.6. 

We therefore decided to match the data to our prediction and extract a common value for 
k. This neglects the possibility that k may show a weak energy dependence over the range 
measured. By combining the results of Fig.||](a) and Fig.|3|, we can plot our prediction 
of cr^ajj. vs. y^i. To obtain k, we leave it as a free parameter and make a y 2 fit to the 
rescaled high energy cosmic ray data for cr^^ r . In Fig.|5] we have replotted the published 
high energy data for cr 1 ™ 1 ^, against our prediction of cr'™^ vs. y/s, using the common value 
of k = 1.33 ± 0.04 ± 0.026, obtained from a x 2 fit- The error of 0.04 is the statistical error of 
the x 2 fit, whereas the error of 0.026 is the systematic error due to the error in the prediction 
of cr^iir- Clearly, we have an excellent fit, with good agreement between AGASA and Fly's 
Eye. The analysis gives x 2 — 1-75 for 6 degrees of freedom (the low x 2 is probably due to 
overestimates of experimental errors). Of course, the improved fit of the cosmic ray data 
has no effect on the fit to the accelerator data. 

Our result for k is interesting — it is close to the value of 1.2 obtained using the SIBYLL 
simulation || for inclusive particle production. This represents a consistency check in the 
sense that our model for forward scattering amplitudes and SYBILL share the same under- 
lying physics. The increase of the total cross section with energy to a black disk of soft 
partons is the shadow of increased particle production which is modeled by the production 
of (mini)-jets in QCD. The difference between the k values of 1.20 and 1.33 could be under- 
stood because the experimental measurement integrates showers in a relatively wide energy 
range, which tends to increase the value of k. 

We predict 108.0 ± 3.4 mb for the total cross section at LHC energy (14 TeV), where 
the error is due to the statistical errors of the fitting parameters. In the near term, we look 
forward to the possibility of repeating this analysis with the higher statistics of the HiRes || 
cosmic ray experiment that is currently in progress and the Auger |10[ Observatory. 
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FIG. 1. B dependence on the pp total cross section a pp . The five curves are lines of constant 
(T p^air) °f 440, 490, 540, 590 and 640 mb — the central value is the published Fly's Eye value, and 
the others are ±1<t and ±2er. The dashed curve is a plot of our QCD-inspired fit of B against <j pp . 
The dot is our value for ^fs = 30 TeV, the Fly's Eye energy. 
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FIG. 2. The simultaneous QCD-inspired fit of total cross section a pp , p and B vs. y/s, in GeV, 
for pp (squares) and pp (circles) accelerator data: (a) a pp , in mb, (b) p, (c) Nuclear slope B, in 
GeV" 2 
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FIG. 3. A plot of the predicted total pp cross section a pp , in mb vs. the measured p-air cross 
section, cr^^, in mb. 



7 



220 
200 
180 
160 
140 
120 
100 
80 
60 
40 
20 




■ cj(pp) 

• o(pbar-p) 

♦ Fly's Eye Air Shower Data 
A AGASA Air Shower Data 

T Aglietta et al, Air Shower Data 



1e+1 



1e+2 1e+3 
Vs, in GeV 



1e+4 



1e+5 



FIG. 4. A plot of the QCD-inspired fit of the total nucleon-nucleon cross section a pp , in mb 
vs. y/s, in Gev. The cosmic ray data that are shown have been converted from c^air to a pp using 
the results of Fig. |3[ 
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FIG. 5. A x 2 fit of the measured AGASA and Fly's Eye data for c^f!^, in mb, as a function 



of the energy, y/s, in GeV. The result of the fit for the parameter k in Eq. (||) is k = 1.33 db 0.04. 
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